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Using  a  first-principles  approach  to  the  ionization  rate,  we  re-examine  some  of  the  prejudices 
concerning  impact  ionization  and  offer  a  new  view  of  the  role  of  thresholds.  We  also  discuss 
trends  of  ionization  coefficients  related  to  the  energy  band  structure  for  silicon  and  a  number  of 
III-V  compounds. 


I.  INTRODUCTION 

Impact  ionization  (the  exact  inverse  of  the  Auger  pro¬ 
cess)  in  semiconductors  has  been  studied  extensively  since 
the  first  theoretical  investigations  of  hot  carrier  transport. 
It  is  now  generally  agreed  that  an  accurate  treatment  of  the 
ionization  process  itself  requires  a  full  band  structure 
model  and  the  associated  complexities.  In  addition,  the 
large  phase  space  relevant  for  the  two-particle  inverse  Au¬ 
ger  process  makes  an  analytic  solution  impossible,  and 
only  detailed  numerical  calculations  can  offer  an  idea  of 
the  magnitude  of  the  matrix  elements  for  the  process.  The 
earliest  attempts  to  understand  impact  ionization  have  by¬ 
passed  these  difficulties  and  instead  focused  on  simple  “in¬ 
tuitive”  concepts  such  as  an  electron  temperature  model, 
the  lucky  electron  models  of  Shockley,'  and  others,^’^  and 
various  models  for  important  parameters  of  the  theory  or 
what  has  been  perceived  to  be  an  important  parameter. 
Prominent  among  these  are  the  phonon  mean-free  path, 
the  ionization  threshold  (the  energy  below  which  impact 
ionization  can  not  occur),  and  the  average  energy  of  net 
emitted  phonons.  In  addition,  one  needs  various  parame¬ 
ters  of  the  energy  band  structure  which,  of  course,  influ¬ 
ence  both  threshold  and  phonon  scattering.  It  was  shown 
subsequently  by  Shichijo  and  Hess^  that  the  energy  depen¬ 
dence  of  phonon  scattering  is  significant  and  its  variations 
with  energy  influence  the  statistics  of  electrons  enor¬ 
mously.  One  therefore  must  include  a  complete  energy 
band  structure  or  at  least  a  strongly  energy-dependent  den¬ 
sity  of  states  to  obtain  the  correct  attempt  frequency  of 
electrons  above  the  threshold  to  impact  ionize.  To  obtain 
rates  of  impact  ionization  one  needs  to  multiply  this  at¬ 
tempt  frequency  by  the  probability  of  the  ionization  event. 
This  probability  was  virtually  exclusively  obtained  from  a 
formula  derived  by  Keldysh.’  This  formula  cannot,  unfor¬ 
tunately,  account  for  the  complexities  introduced  by  the 
Bloch  functions  or  an  energy  band  structure  with  its  min¬ 
ima  at  more  than  one  symmetry  point  (typically  at  T,  X, 
and  L  in  the  technologically  important  semiconductors). 
The  problems  of  the  Keldysh  formula  were  first  signaled  by 
a  numerical  calculation  of  Kane^  which  showed  that  the 
ionization  rate  of  silicon  increases  above  the  threshold 
much  more  softly  than  predicted  by  Keldysh.  Another  set 
of  difficulties  becomes  apparent  from  discussions  by  Ca- 
passo  regarding  the  uncertainty  principle.^  Since  the  elec¬ 
tron  system  is  in  contact  with  a  phonon  bath,  it  does  not 
form  a  Hamiltonian  system,  and  therefore  does  not  have  a 


strict  threshold  at  all.  This  was  quantitatively  shown  by 
Bude  et  al}  who  also  showed  that  the  high  electric  fields 
which  are  typically  present  when  ionization  becomes  im¬ 
portant  have  a  significant  influence  on  threshold  and  the 
softness  of  the  rate.  We  show  that  as  a  consequence  of 
these  facts,  many  of  the  theoretical  discussions  in  journals 
and  textbooks’'"  are  misleading  or  even  completely 
wrong. A  typical  example  is  the  calculation  of  threshold 
energy  using  the  effective  masses  of  the  band  edges. 

Questions  such  as  the  dependence  of  ionization  rates 
on  the  crystallographic  direction  need  even  more  careful 
consideration.  As  an  example,  we  will  discuss  the  k- 
dependent  ionization  thresholds,  revisiting  the  treatment  of 
Anderson  and  Cmwell.'^ 

Our  findings  show  that  the  concept  of  hard  thresholds 
for  the  impact  ionization  rate  needs  to  be  abandoned  for  all 
materials  which  we  have  investigated.  Our  numerical  sim¬ 
ulations  have  instead  shown  effective  thresholds  which  are 
significantly  higher  in  the  bands  than  the  laws  of  energy 
and  momentum  conservation  dictate.  The  existence  of  a 
threshold  by  itself  cannot  be  deduced  with  physical  rigor 
since  the  system  of  electrons  is  not  Hamiltonian.  The  pho¬ 
non  bath  plays  an  important  role  and  broadens  the  thresh¬ 
old.*  Also,  strong  electric  fields,  usually  present,  exert  a 
similar  influence  due  to  the  intracollisional  field  effect. 

An  even  more  important  factor  for  the  effective  thresh¬ 
old  is  the  functional  form  of  the  density  of  states  in  the 
respective  semiconductor.  The  minima  of  highest  density 
of  states  (usually  at  20  play,  therefore,  a  significant  role 
for  the  location  of  the  effective  threshold.  The  density  of 
states  also  deeply  influences  the  magnitude  of  the  average 
energy  of  the  electrons.  This  magnitude  in  turn  determines 
whether  the  ionization  occurs  due  to  electrons  of  the 
“bulk”  or  “center”  of  the  distribution  function  or  due  to  a 
few  lucky  ones  as  postulated  by  Shockley.  The  correlation 
of  effective  threshold  and  average  energy  to  the  density  of 
states  works  against  Shockley’s  lucky  electron  theory  as  far 
as  the  impact  ionization  coefficients  (a  and  /?)  are  con¬ 
cerned.  The  functional  dependence  of  a  and  p  on  the  field 
{exp(c//')  or  exp(c//'^)]  is  also  closely  correlated  to  the 
rise  of  the  density  of  states  in  the  important  energy  range 
and  not  as  much,  as  is  claimed  in  most  textbooks  (includ¬ 
ing  Ref.  14),  on  the  question  of  whether  lucky  electrons 
are  important  for  the  ionization  or  not.  We  realize  that  our 
findings  represent  a  very  significant  deviation  from  the  con¬ 
ventional  picture  of  impact  ionization  in  semiconductors. 
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Our  claims  are  not  absolute  in  the  sense  that  we  can  back 
them  by  mathematical  proof.  We  can  prove  them  only  by 
computer  experiments  for  a  limited  but  important  class  of 
semiconductor  materials.  We  believe  that  our  computer 
experiments,  which  combine  and  optimize  the  numerical 
approaches  of  Kane,*’  BarafF,'’  and  Hess  and  Shichijo,*  are 
the  most  complete  hitherto  preformed.  They  still  are  not 
based  on  exact  band  structures  or  scattering  rates  but  on 
what  one  obtains  within  the  framework  of  empirical 
pseudopotentials.  This  fact,  however,  should  not  influence 
the  “bigger  picture”  that  is  discussed  here. 

II.  THE  WAVE  VECTOR  DEPENDENT  THRESHOLD 
ENERGY  AND  RESULTS  FOR  THE  RATE  IN 
SILICON 

The  first  paper  to  carefully  consider  extending  the 
threshold  energy  concept  to  realistic  band  structures  sug¬ 
gested  that  because  of  the  anisotropy  and  the  multiplicity 
of  bands  in  real  crystals,  the  threshold  energy  should  be 
defined  as  a  function  of  direction  in  the  crystal.'^  An  elec¬ 
tron  with  a  given  initial  wave  vector,  k,,  could  impact 
ionize  only  if  i^  energy  exceeded  a  minimum  energy, 
J?th(^i),  where  A'l  is  the  unit  normal  in  the  kj  direction. 
Here,  we  review  the  approach  used  by  Anderson  and 
Crowell’’  to  calculate  £,(,(k|),  and  demonstrate  an  incon¬ 
sistency  which  exaggerates  the  importance  of  wave  vector 
dependent  thresholds  of  the  wave  vector  anisotropy  of  the 
ionization  rate. 

First,  we  define  the  particles  involved  in  the  electron 
ionization  process;  considerations  for  the  hole  process  fol¬ 
low  in  the  same  manner.  Let  k]  be  the  initial  wave  vectors 
for  a  conduction  band  electron  which  collides  with  a  va¬ 
lence  band  electron  of  wave  vector  k2.  The  two  final  states 
are  conduction  band  states  of  wave  vector  k3  and  k4.  In  the 
semiclassical  approximation,  energy  is  conserved.  Further¬ 
more,  the  matrix  element  for  the  electron-electron  interac¬ 
tion  enforces  crystal  momentum  conservation,  and  we  can 
write 

k|=k3-|-k4  —  k2-l-G,  (1) 

£„,(k,)=£'„3(k3)+£„^(k4)-£'„^(k2),  (2) 

where  G  is  a  reciprocal  lattice  vector  and  £„(k)  is  the 
energy  of  an  electron  in  band  n  with  wave  vector  k. 

Given  an  initial  k],  one  can  search  ail  possible  sets  of 
k3,  k4,  k2,  and  G  satisfying  ( 1 )  and  ( 2 )  to  find  the  mini¬ 
mum  energy  that  the  incident  particle  could  have,  E„(V.^). 
If  the  actual  energy  of  the  initial  state  £„^(ki)  is  greater 
than  £„,(ki),  then  the  impact  ionization  rate  could  be  fi¬ 
nite  for  this  carrier.  If  £„,(ki)  >  £„(ki)  for  all  n,,  then 
impact  ionization  is  impossible  for  a  carrier  with  wave  vec¬ 
tor  k).  is  defined  as  follows;  for  a  given  direction  in 

the  Brillouin  zone,  k,  E„,  and  E„^  are  compared.  The  pos¬ 
sible  values  of  E^^(k)  are  given  by  the  energies  where 
E„{k)  =  £„|(k)  as  k  varies  along  the  direction  fc.”  Then, 

E,h(*)smin(£'„(k);£’,„(k)=£'„|(k), 


£„(k)=  min  [E^^lkj) -|-£„^(k4) -£.  (k2), 

C.k^yf^y  =  2,4 

k3-|-k4-k2-hG  =  k].  (3) 

An  application  of  the  Lagrange  multiplier  rule  to  the 
constrained  minimum  in  (3)  requires  that  the  group  ve¬ 
locities  of  states  2,3,  and  4  be  equal. 

V£„^(k2)=V£„^(k3)=V£„^(k4).  (4) 

This  is  a  necessary  condition  for  minimization.  A  mini¬ 
mum  is  guaranteed  to  exist  since  E„  and  £„^  are  bounded 
from  below  and  £„,  is  bounded  from  above.  Of  all  the 
states  2,  3,  and  4  for  which  the  group  velocities  are  equal, 
it  is  simple  to  choose  the  one  which  corresponds  to  the 
absolute  minimum  required  in  (3). 

Although  the  above  definition  of  £th(^)  describes 
whether  or  not  an  electron  in  a  certain  direction  can  im¬ 
pact  ionize,  the  behavior  of  the  rate  beyond  this  threshold 
is  unclear.  If  one  assumes  that  it  increases  rapidly  beyond 
the  threshold,  then  this  theory  provides  an  extension  of  the 
energy-dependent  thresholds  of  Keldysh  and  Shockley. 

The  problem  as  stated  above  has  not  been  solved  cor¬ 
rectly  in  the  original  paper  by  Anderson  and  Crowell.” 
They  propose  that  in  evaluating  £„  for  a  given  k,  the  min¬ 
imization  over  k2,  k3,  and  k4  is  achieved  by  states  which  are 
parallel: 

k2  k3  k*  /s 

Consequently,  their  search  for  solutions  of  (4)  extends 
only  over  the  subset  of  states  2,  3,  and  4  satisfying  (5). 
Since  their  search  range  is  drastically  reduced,  solutions 
for  £,h(k)  were  easily  accomplished  numerically.  The  ac¬ 
tual  solution  of  these  equations  requires  considerably  more 
effort. 

Clearly,  the  solution  for  E„  requires  a  minimization 
over  all  values  of  k2,  k3,  and  k4  for  a  given  k.  Only  if  the 
bands  are  isotropic  is  the  solution  for  £„,  achieved  by  par¬ 
allel  states.  To  demonstrate  this  we  propose  the  following 
simple  example.  Instead  of  using  spherical-parabolic 
bands,  we  take  ellipsoidal-parabolic  bands  in  two  dimen¬ 
sions; 


Ec=a^  kl+Uyk^, 

(6) 

E„=~Eg-P(kl+kl), 

(7) 

where  E^  and  £„  correspond  to  the  conduction  and  valence 
bands,  respectively,  and  Eg  to  the  band  gap  energy.  Apply¬ 
ing  (4)  and  kj-|-k4  — k2  =  k  to  the  bands  defined  in  (6) 
requires  that  the  minimum  in  £„(k)  be  achieved  by  states 
satisfying 


^  =  ^=15(0;, -I- 2/3)  'k^. 

(8) 

=P(ay+W)  ~'ky, 

(9) 

k,=[2P{a,  +  W)-'-\]k,, 

(10) 
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TABLE  1.  Wave  vector  dependent  threshold  in  silicon. 


Thres- 

Direc- 

hold 

lion 

(eV) 

ki 

k4 

Anderson  and  Crowell 


<100) 

I.l 

-0.40(1,0,0) 

0.00(1,0,0) 

0.80(1,0,0) 

0.80(1,0,0) 

(110) 

2.1 

-0.47(1,1,0) 

0.19(1,0,0) 

0.86(1,1,0) 

0.86(1,1,0) 

(111) 

3.1 

-0.02(1,1,1) 

0.00(1,1,1) 

0.99(1,1  1) 

0.99(1,1,1) 

Czajkowski  et 

al 

(100) 

1.18 

-0.36(1,0,0) 

0.01(1,0,0) 

0.82(1,0,0) 

0.82(1,0,0) 

(110) 

2.10 

-0.69(1,1,0) 

0.31(1,0,0) 

1.22(1,1,0) 

1.42(1,1,0) 

(111) 

3.39 

-0.07(1,1,1) 

0.01(1,1,1) 

0.84(1,1,1) 

0.84(1.1,1) 

^7. 

=  (2^( 

ay+2E)~'- 

-I]*,- 

(11) 

These  minimizing  states  are  parallel  to  the  initial  state  k 
only  when  Therefore,  in  general,  for  complicated 

band  structures  there  is  no  simple  relationship  between  the 
directions  of  the  minimizing  states  and  the  initial  state. 

In  silicon,  in  which  the  minimum  of  the  conduction 
bands  is  at  the  X  point,  this  error  is  particularly  bad.  Listed 
in  Table  I  art  typical  results  (Anderson  and  Crowell'^  and 
Czajkowski  et  al.'^)  for  the  lowest  thresholds  in  several 
symmetry  directions  in  Si  with  wave  vectors  of  the  “thresh¬ 
old”  states. 

According  to  Table  I,  the  threshold  in  the  (Ill)  di¬ 
rection  in  Si  is  greater  than  3.0  eV,  and  therefore,  the  only 
initial  states  in  the  ( 1 1 1 )  direction  which  can  impact  ionize 
must  have  an  energy  of  at  least  3  eV.  In  Fig.  1  we  show  our 
results  which  demonstrate  that  it  is  clearly  possible  for 
electrons  along  the  <  1 1 1  >  direction  with  energies  much  less 
than  3  eV  to  impact  ionize.  Figure  2  shows  the  pseudopo¬ 
tential  band  structure  in  these  crystallographic  directions. 
The  ionization  rate  is  nonzero  and  large  in  the  (111)  di¬ 
rection  even  at  electron  energies  as  low  as  2.2  eV.  This 
directly  contradicts  the  results  shown  in  Table  I  which 


FIG.  I .  Ionization  rates  in  silicon  for  the  two  lowest  conduction  bands  as 
a  function  of  wave  vector.  Upper  curves,  second  band;  lower  curves,  first 
band.  These  results  are  obtained  as  described  in  Ref  8  and  are  in  agree¬ 
ment  with  Kane’s  theory  (Ref  6). 


Wav»-vaetor  (2ic/a| 


FIG.  2.  Band  structure  of  silicon  for  two  lowest  conduction  bands  as  a 
function  of  wave  vector.  Upper  curves,  second  band;  lower  curves,  first 
band. 


suggest  that  in  the  (111)  direction,  only  electrons  whose 
energy  is  greater  than  3.1  eV  can  impact  ionize.  The  same 
is  true  of  the  ionization  rate  in  the  (110)  direction.  The 
thresholds  derived  in  the  simplified  manner  of  Anderson 
and  Crowell  are  therefore  inaccurate.  Note  that  in  all  sub¬ 
sequent  publications  known  to  the  authors,  it  was  this  sim¬ 
plified  inaccurate  model  of  thresholds  that  has  been  used. 
In  addition,  the  wave  vector  dependent  thresholds  in  Eq. 
(3)  can  tell  little  about  the  actual  behavior  of  the  ioniza¬ 
tion  rate  for  energies  above  £,|,  since  they  distinguish  only 
if  a  particular  electron  can  ionize  as  a  function  of  its  direc¬ 
tion  and  tell  nothing  about  the  actual  behavior  of  the  ion¬ 
ization  rate  for  energies  above  E^^(k). 

We  would  also  like  to  reiterate  that  a  quantum  treat¬ 
ment  of  the  electrons  and  the  phonon  bath  does  not  result 
in  any  absolute  threshold  at  all  due  to  collision  broadening 
and  the  intracollisional  field  effect.  In  the  following,  we 
ignore  this  “smearing  out”  and  show  that  the  semiclassical 
inclusion  of  band  structure  and  phonon  scattering  alone  is 
sufficient  to  justify  the  concept  of  an  effective  threshold 
which  is  typically  even  higher  than  that  of  Anderson- 
Crowell  and  which  shows,  in  the  cases  which  we  have 
considered  in  our  computer  experiments,  very  little  anisot¬ 
ropy  with  respect  to  the  crystallographic  direction. 


ill.  IMPACT  IONIZATION  IN  POLAR 
SEMICONDUCTORS 

In  this  section,  we  examine  the  average  energy- 
dependent  ionization  rates  of  several  important  III-V  com¬ 
pounds  and  the  influence  of  these  rates  on  bulk  transport. 
We  use  the  same  theory  as  developed  for  silicon.^’*  Of 
special  importance  for  our  considerations  is  the  average 
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ionization  rate  for  electrons  of  a  given  energy  which  can  be 
expressed  as 

2Jgzd^k8[E„(k)-E] 


lj8zd^kR„(k)S[E„(k)-E] 

G(E) 


(12) 


where  R„(k)  is  the  ionization  rate  (transition  probability 
per  unit  time)  for  an  electron  with  wave  vector  k  and  band 
index  n.  The  total  density  of  states  at  energy  E  is  G{E), 
and  f  B2  d^k  is  an  integration  over  the  first  Brillouin  zone. 
Our  theory  for  R„(k)  is,  as  mentioned,  very  similar  to  the 
theory  of  Kane,  and  has  been  described  in  Refs.  6  and  8 
extensively.  Below,  we  review  some  of  the  salient  features 
for  this  theory  and  present  the  results  of  our  computer 
experiments. 

A.  Energy  dependent  ionization  rates 

Two  quantities  strongly  influence  the  energy- 
dependent  ionization  rate:  the  threshold  energy  itself,  and 
its  hardness  or  softness  (how  impact  ionization  and  pho¬ 
non  scattering  compete).  The  Keldysh  formula,^ 

P  /£-£:, h\^ 

(  £.h  )’ 

(13) 

has  two  adjustable  parameters  for  a  given  phonon  scatter¬ 
ing  rate  Tp^,  — £th,  which  is  the  threshold  energy,  and  P 
which  is  a  measure  of  the  hardness  of  the  threshold.  (Ac¬ 
tually,  Pis  probably  a  better  measure  of  the  hardness  of  the 
threshold  if  the  phonon  energy  relaxation  rate  is  used  in 
place  of  Tph.'^)  Monte  Carlo  (MC)  calculations,  which  use 
the  Keldysh  formula,  have  deduced  a  relatively  hard  ion¬ 
ization  threshold  for  most  III-V  semiconductors  such  as 
GaAs,  and  a  soft  threshold  for  Si.^  '*"^'  Although  used  in 
most  high  field  transport  simulations,  the  Keldysh  formula 
applies  only  to  materials  with  parabolic  bands  and  a  direct 
band  gap.  Real  materials  deviate  significantly  from  this 
simple  model  band  structure.  It  also  does  not  represent  a 
good  parametrization  of  a  general  ionization  rate  since  it 
allows  only  for  a  quadratic  dependence  on  energy. 

A  direct  calculation  of  (12)  entails  considerable  com¬ 
putational  effort,  and  has  until  now  been  carried  out  only 
for  Si.*  Here,  we  present  for  the  first  time  the  energy  de¬ 
pendent  ionization  rates  for  the  HI-V  materials  GaAs, 
InAs,  InP,  and  the  alloy  Gao ^alno.jvAs.  The  calculations 
have  used  a  model  similar  to  the  Kane  model  for  impact 
ionization*  for  the  semiclassical  ionization  rates  in  which 
local  empirical  pseudopotential  methods  were  used  to  cal¬ 
culate  the  band  structure  and  wave  functions  for  each  ma¬ 
terial. Band  structure  fits  were  made  to  experimentally 
determined  quantities  as  referenced  in  Refs.  21  and  24. 
Representative  band  structure  parameters  for  these  mate¬ 
rials  are  shown  in  Table  II.  The  density  of  states  (DOS) 
for  each  material  is  plotted  in  Fig.  3. 

There  is  one  inconsistency  in  our  Gao43lno  57AS  band 
structure.  While  most  parameters  such  as  the  band  gap  and 
the  elFective  masses  are  easily  fit,  the  T  to  L  intervalley 


Table  II.  Band  structure  parameters. 


Material 

Er-i 

Er-x 

mf 

GaAs 

1.4 

0.25 

0.32 

0.067 

InAs 

0.32 

1.10 

1.75 

0.032 

InP 

1.3 

0.75 

1.20 

0.075 

57  As 

0.75 

0.78 

1.10 

0.041 

separation  turns  out  to  be  about  0.78  eV,  a  typical 
result  for  pseudopotential  and  kp  theory,^'  whereas  some 
measurements  have  shown  a  separation  of  0.55  eV.^’  In 
general,  the  pseudopotential  minimum  separation  is  not  as 
precise  as  desirable  (see  the  results  for  GaAs  in  Table  II). 
Also,  the  split-off  bands  have  not  been  accurately  modeled 
using  local  pseudopotentials.  However,  they  contribute 
very  little  to  the  ionization  rate  because  of  their  relatively 
small  density  of  states.  Therefore,  the  local  pseudopoten¬ 
tials  without  spin-orbit  effects  should  give  a  fairly  accurate 
estimation  of  the  ionization  rate,  particularly,  if  one  keeps 
in  mind  that  the  valley  separation  energy  should  be  some¬ 
what  shifted. 

In  going  from  GaAs  to  InAs,  the  band  gap  shrinks 
from  1.4  to  0.32  eV  while  the  T-X  and  T-L  intervalley 
separations  rise  dramatically.  This  leads  to  a  very  small 
density  of  states  for  low  energies  in  InAs  and 
Gao4jIno57As  compared  to  that  in  GaAs.  However,  InP 
lies  between  these  extremes  in  that  it  has  a  large  band  gap 
but  the  P-T’  and  P-L  intervalley  separations  are  rather 
large.  These  band  structure  features  lead  to  two  very  dis¬ 
tinct  behaviors  for  the  bulk  ionization  coefficient  as  will  be 
shown  below. 

Shown  in  Fig.  4  are  the  calculated  results  for  the 
energy -dependent  ionization  rates  for  GaAs,  InAs,  InP, 
and  Gao  43lno  57AS.  Plotted  along  with  each  rate  for  com¬ 
parison  are  Keldysh  formula  ionization  rates,  which  are 
typical  of  those  used  in  Monte  Carlo  simulations”  (param¬ 
eters  taken  from  Ref.  21 ). 


1*10“ 


FIG.  3.  DOS'  solid  line,  GaAs:  dashed  line,  InAs;  dot-dashed  line:  InP, 
and  dotted  line:  GaojjlnojiAs. 
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FIG.  4,  Ionization  rates:  upper  set  of  three  curves.  Keldysh  formula: 
lower  set  of  four  curves,  calculated  ionization  rates.  Solid  line;  GaAs, 
dashed  line:  InAs,  dot-dashed  line:  InP,  and  dotted  line:  Gao  43lno  57 As. 


The  semiclassical  ionization  rate  thresholds  for  these 
materials  are  always  larger  than  the  band  gap  energy.  In 
GaAs  for  instance,  the  ionization  rate  becomes  nonzero  for 
energies  about  0.5  eV  larger  than  the  band  gap,  wherejis  in 
the  InAs-like  compounds,  the  threshold  is  closer  to  the 
band  gap  energy  (0.1  eV  larger  than  the  band  gap  energy). 
The  threshold  m  InP  is  almost  1  eV  larger  than  the  band 
gap  energy.  However,  the  actual  energy  at  which  the  ion¬ 
ization  rate  becomes  nonzero  is  difficult  to  determine  when 
using  a  Monte  Carlo  integration  algorithm  as  in  Ref.  6. 
This  mathematical  threshold  is  also  largely  irrelevant  as 
discussed  above.  A  more  practical  measurement  of  an  ion¬ 
ization  threshold  is  to  define  the  threshold  as  the  minimum 
energy  an  electron  must  attain  to  contribute  a  large  per¬ 
centage  to  a  macroscopic  quantity  such  as  the  impact  ion¬ 
ization  coefficient  a.  This  will  be  defined  more  clearly  be¬ 
low  in  the  discussion  on  bulk  transport  properties. 

For  realistic  band  structures,  there  is  a  complicated 
interplay  between  the  threshold  energy  and  hardness  of  the 
threshold.  It  is  possible,  however,  to  make  the  following 
rough  classifications  pertaining  to  direct  gap  materials.  For 
small  band  gap  materials  with  large  Ej-.x  and  Ej-.l,  the 
threshold  energy  is  close  to  the  band  gap  energy  and  the 
threshold  is  soft.  In  this  case,  electrons  with  an  energy 
slightly  greater  than  the  band  gap  energy  above  the  con¬ 
duction  band  minimum  are  located  in  the  F  minimum.  The 
final  states  are  also  located  in  the  F  minimum  where  the 
density  of  states  is  very  small.  In  these  cases,  the  ionization 
rate  will  be  roughly  proportional  to  the  product 
G^(Ei)G^(E4),  where  Ej  and  E^  are  the  final  state  energies. 
Therefore,  the  ionization  rate  near  threshold  will  increase 
slowly  until  it  is  possible  for  the  intervalley  states  to  be¬ 
come  final  states.  We  have  numerically  shown  that  this  is 
the  case  for  InAs  and  Gag^jlno  57AS. 

For  larger  band  gap  materials  with  small  E^.x  and 
Ef^.i,  the  threshold  energy  will  be  slightly  larger  than  the 
band  gap  and  the  threshold  will  be  harder.  If  E^-.x  and  E^.i^ 
are  small,  the  low  energy  DOS  will  be  large  and  will  in- 
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crease  rapidly  with  energy;  hence,  the  available  final  states 
for  impact  ionization  increase  rapidly  with  initial  electron 
energy.  Although  electrons  near  threshold  will  lie  high  in 
the  band  structure  and  will  have  a  wide  range  of  wave 
vector  magnitudes,  the  low-energy  final  states  have  a  large 
wave  vector  spread,  and  energy  and  momentum  conserva¬ 
tion  raises  the  threshold  energy  only  a  little  above  the  band 
gap  energy.  This  is  true  for  GaAs. 

If  the  band  gap  is  large  and  Ej^  x  z.  a*"®  large, 

the  energy  and  momentum  conservation  requirement  can 
force  the  threshold  energy  to  be  substantially  larger  than 
the  band  gap  energy.  Since  the  final  states  are  restricted  to 
the  F  valley  until  high  initial  particle  energies  are  reached, 
it  may  be  impossible  to  conserve  both  energy  and  momen¬ 
tum  unless  the  initial  state  electron  energy  is  much  larger 
than  the  band  gap  energy.  This  condition  pertains  to  InP 
with  a  1.3  eV  band  gap  and  a  valley  separation  of  more 
than  0.7  eV. 

In  addition,  the  Coulomb  matrix  element  responsible 
for  the  electron-electron  interaction  can  play  an  impiortant 
role  in  determining  the  ionization  rate  in  III-V’s,  especially 
for  energies  near  threshold.  Since  the  matrix  element  has  a 
rough  l/q^  dependence  (see  Refs.  6  and  8),  the  ionization 
rate  for  initial  state  energies  for  which  the  average  g  ex¬ 
changed  is  small  are  enhanced  by  factors  of  10  to  100 
compared  to  ones  for  which  9  is  a  substantial  fraction  of 
the  BZ.  Because  the  average  q  exchanged  is  much  smaller 
for  energies  near  threshold,  the  average  matrix  elements 
for  these  states  are  larger  than  those  of  higher  energy.  This 
effect  is  most  pronounced  for  smaller  band  gap  materials 
such  as  InAs. 


B.  Bulk  transport  and  ionization  coefficient  a 

In  all  four  materials  which  we  have  considered  here, 
the  ionization  raves  calculated  are  much  smaller  than  those 
obtained  by  the  empirical  Keldysh  formula  which  has  been 
used  in  past  Monte  Carlo  simulations  (see  Fig.  4).  How¬ 
ever,  it  is  possible  to  obtain  Monte  Carlo  fits  using  the  new 
Kane-type  results  to  experimentally  obtained  values  of  the 
impact  ionization  coefficient  a  without  appreciably  chang¬ 
ing  the  total  phonon  scattering  rate.  We  have  approxi¬ 
mated  the  impact  ionization  coefficients  a  for  electrons  and 
j3  for  holes  as 

2„/Bzrf'fc[l-A,(k)]Ravel£„„(k)] 

^nSBzdM\-f„W]v,  ’ 

where  u„(k)  is  the  group  velocity,  f„(k)  is  the  distribution 
function,  and  and  are  the  conduction  and  valence 
band  indices,  respectively.  The  fact  that  one  can  virtually 
obtain  the  same  result  for  a  using  two  totally  different  sets 
of  ionization  rates  has  been  noticed  by  Tang  for  the  case  of 
silicon,’*  and  represents  a  general  property  of  the  Boltz¬ 
mann  equation  and  scattering  mechanisms  that  involve 
some  form  of  threshold  (see,  for  instance  Ref.  26). 
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FIG.  5.  Phonon  scattering  rates  used  in  MC  simulations.  Solid  line: 
GaAs,  dashed  line:  InAs,  dot-dashed  line:  InP,  and  dotted  line: 

GA)43^tto57As. 


We  have  performed  full  band  structure  bulk  Monte 
Carlo  simulations  and  have  calculated  distribution  func¬ 
tions  and  macroscopic  quantities  such  as  a  and  drift  veloc¬ 
ity  as  a  function  of  field  for  four  III-V  compounds.  Our 
Monte  Carlo  simulation  uses  two  pseudopotential  conduc¬ 
tion  hands  for  drift  due  to  the  field,  and  includes  several 
scattering  mechanisms  such  as  polar-optical  phonon  scat¬ 
tering  and  intervalley  scattering  (see  Ref.  14).  The  optical 
phonon  energies  we  have  used  are  similar  to  those  of  Refs. 
19-21.  In  addition,  the  impact  ionization  rates  shown  here 
were  also  included  as  energy-dependent  scattering  mecha¬ 
nisms.  The  phonon  scattering  rates  at  high  fields  are  taken 
to  be  proportional  to  the  total  density  of  states  in  the  con¬ 
duction  band,  which  is  a  good  approximation  for  deforma¬ 
tion  potential  scattering.'® 

Figure  5  shows  the  total  phonon  scattering  rates  used 
in  these  MC  simulations.  These  scattering  rates  where  ad¬ 
justed  to  roughly  match  bulk  measurements  of  a  for 
GaAs,^^  InP,^®’^’  and  Gao.43lno  57As,®°  but  because  there  is 
little  data  available  on  the  InAs  ionization  coefficient,  pho¬ 
non  scattering  rates  comparable  to  Gao43lno57As  were 
used  in  this  case.  Figure  6  shows  the  values  of  a  corre¬ 
sponding  to  these  phonon  scattering  rates. 

Although  bulk  transport  parameters  can  be  roughly 
matched  using  both  the  Keldysh-type  and  our  new  Kane- 
type  results,  the  distribution  functions  which  result  are 
quite  different  and  therefore,  a  different  interpretation  of 
the  ionization  threshold  is  in  order.  For  instance,  in  GaAs, 
typically  used  phonon  coupling  constants  predict  total 
phonon  scattering  rates  which  peak  at  3x  lO'^  s~’  near  1.5 
and  3.3  eV,  have  an  average  high  energy  value  of  about 
2x  lO'^  s“'  (see  Ref  21 ),  and  an  ionization  rate  from  the 
Keldysh  formula  as  shown  in  Fig.  4.  In  these  models,  the 
ionization  rate  quickly  surpasses  the  total  phonon  scatter¬ 
ing  rate  near  2  eV.  Because  both  the  energy  lost  in  impact 
ionization  and  the  ionization  rate  are  large,  almost  no  elec¬ 
trons  can  achieve  energies  greater  than  2.25-2.5  eV.  There¬ 
fore,  in  these  GaAs  simulations,  the  distribution  functions 


FIG.  6.  Calculated  impact  ionization  coefficient  a  as  a  function  of  inverse 
field,  \/F.  Solid  line:  GaAs,  dashed  line:  InAs,  dot-dashed  line.  InP,  and 
dotted  line:  Gao43lnfl57As. 


have  truncated  high-energy  tails,  and  a  is  determined  pri¬ 
marily  by  how  many  electrons  can  reach  an  effective 
threshold  (near  2.0  eV).  This  number  is  critically  deter¬ 
mined  by  the  phonon  scattering  rate  and  qualitatively  be¬ 
haves  corresponding  to  the  Shockley  "lucky  electron” 
model.  The  same  behavior  is  seen  in  other  III-V  materials 
when  simulated  in  this  manner.'’  Notice,  however,  that  in 
order  to  fit  experiment,  the  Keldysh  parameters  for  each 
materia]  vary  by  orders  of  magnitude  from  material  to 
material  (see  Fig.  4),  and  phonon  rates  can  vary  by  factors 
of  3  to  4. 

In  contrast,  a  different  picture  of  ionization  thresholds 
is  obtained  from  the  ionization  rates  calculated  here  in  the 
Kane-type  approach.  For  example,  in  GaAs,  the  total  pho¬ 
non  scattering  rate  is  roughly  20%  lower  than  that  of  Ref 
21,  but  the  ionization  rates  are  a  factor  of  1000  lower.  In 
our  simulations,  the  ionization  rate  plays  little  role  in  de¬ 
termining  the  distribution  function  up  to  very  high  fields, 
since  it  is  several  orders  of  magnitude  smaller  than  the 
phonon  scattering  rate.  The  ionization  coefficient  a  is  in¬ 
stead  determined  by  the  high-energy  tail  of  the  distribution 
function  above  2.0  eV  integrated  with  the  ionization  rate. 
In  this  model,  the  ionization  rate  samples  the  high-energy 
tail  without  disturbing  it  very  much.  Figures  7  and  8  show 
the  distribution  function  in  GaAs  and  Gao43lno  57As  as  a 
function  of  energy  for  various  values  of  electric  field.  In 
GaAs,  as  the  field  is  increased  and  the  distribution  heats 
up,  however,  there  are  electrons  which  reach  3.5  eV  where 
the  ionization  rate  is  between  lO'^  and  lO'®  s“'.  At  these 
energies,  the  energy  loss  rate  due  to  impact  ionization  ex¬ 
ceeds  that  due  to  phonons,  and  the  part  of  the  high-energy 
tail  above  these  energies  is  truncated  by  impact  ionization. 
However,  impact  ionization  has  little  effect  on  the  distri¬ 
bution  function  for  energies  less  than  the  extremely  high 
value  of  3.5  eV.  The  same  picture  of  impact  ionization 
holds  for  the  other  materials  which  we  have  investigated 
here.  Notice  that  our  results  contradict  another  common 
prejudice  in  the  theory  of  impact  ionization.  It  often  has 
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FIG.  7.  Normalized  electron  distribution  function  in  GaAs  including 
DOS. 


FIG.  9.  Average  electron  energy  as  a  function  of  inverse  field.  Solid  line-. 
GaAs,  dashed  line:  InAs,  dot-dashed  line:  InP.  and  dotted  line: 


been  assumed  that  the  ionization  rate  cannot  be  much 
smaller  than  the  phonon  scattering  rate  to  give  significant 
values  of  a.  Moll^'  considered  the  effect  of  several  phonon 
collisions  in  between  ionizing  events,  but  his  arguments  are 
based  on  assumptions  which  do  not  account  for  the  com¬ 
plicated  form  of  the  energy  distribution  of  the  electrons. 

Soft  thresholds  have  also  been  recently  included  in  the 
“Ridley  lucky  electron  model.”  These  models  incorporate 
Keldysh-type  ionization  rates  whose  magnitudes  are  much 
closer  to  those  calculated  here.'^'^^  These  models  can  suc¬ 
cessfully  predict  many  features  pertaining  to  the  ionization 
coefficients  a  and  p-,  however,  because  they  are  phenome¬ 
nological  models  they  do  not  give  an  accurate  representa¬ 
tion  of  the  ionization  process  and  thus  cannot  predict  some 
important  aspects  of  the  distribution  function. 

A  more  practical  definition  of  thresholds  pertaining  to 
these  very  low  ionization  rates  is  to  define  an  effective 
threshold  energy.  If  we  define  a  “significant”  ionization 
rate  as  a  >  1  and  look  to  see  which  portions  of  the  ioniza¬ 
tion  rate  versus  energy  curve  contribute  to  a>  i,  we  find 


Energy  MMtwed  from  Conduction  Band  Edge  loV) 


FIG.  8.  Normalized  electron  distribution  function  in  Gao43lno,7As  in¬ 
cluding  DOS. 


that  for  these  materials  there  is  no  “significant”  contribu¬ 
tion  to  a  for  values  of  the  ionization  rate  less  than  lO**  s  '. 
The  eflFective  ionization  threshold  is  then  the  lowest  energy 
for  which  the  ionization  rate  is  greater  than  lO’  s~'. 

As  a  result  of  the  small  ionization  rates,  the  average 
energy  of  the  distribution  function  is  much  higher  in  our 
simulation.  Figure  9  shows  the  average  energies  for  the 
four  materials  calculated  here.  These  plots  demonstrate  an 
important  trend.  The  higher  the  intervalley  energies,  the 
higher  the  average  energy.  For  example,  InAs,  which  has 
the  highest  intervalley  energy  separations  exhibits  the  high¬ 
est  average  energy,  whereas  GaAs,  which  has  the  lowest 
intervalley  energy  separations  exhibits  the  lowest  average 
energies.  This  situation  is  expected  to  occur  because  there 
is  little  energy  loss  due  to  impact  ionization  compared  to 
phonon  scattering.  It  is  this  fact  which  leads  to  two  very 
dissimilar  behaviors  in  a  vs  1  plots.  In  GaAs  and  InP,  a 
behaves  as  where  Fq  is  some  large  constant.  For 

these  materials,  a  drops  quickly  to  zero  for  fields  less  than 
200  kV/cm.  This  is  a  result  of  high  effective  threshold 
energies,  particularly  in  InP,  and  low  average  energies  in 
GaAs.  In  contrast,  in  InAs  and  Gao.43lno  57AS,  a  deviates 
significantly  from  the  form  due  to  relatively  low 

threshold  energies  and  high  average  energies. 

IV.  CONCLUSIONS 

We  have  reviewed  the  role  of  thresholds  in  impact  ion¬ 
ization  considering  the  k  dependence  and  energy  depen¬ 
dence  of  the  ionization  cross  section  as  well  as  the  notion  of 
“effective”  ionization  thresholds  for  bulk  transport.  In  our 
discussion  of  the  k  anisotropy,  we  have  shown  that  earlier 
theories  of  ionization  rate  anisotropy  have  been  calculated 
in  error,  and  we  have  presented  a  more  rigorous  calcula¬ 
tion  of  the  ionization  rate  in  silicon  as  a  function  of  k.  This 
more  rigorous  result  demonstrates  that  these  earlier  theo¬ 
ries  have  overestimated  the  anisotropy  of  the  ionization 
rate.  We  have  also  calculated  from  first  principles  the  ion¬ 
ization  rates  in  GaAs,  InAs,  InP,  and  Gao43lno  57AS  which 
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are  found  to  be  much  smaller  than  the  phonon  scattering 
rates  and  exhibit  soft  threshold  behavior.  These  rates  differ 
markedly  from  the  established  Keldysh  formula  ionization 
rates.  The  differences  imply  much  “hotter”  electron  energy 
distribution  functions  and  much  higher  “effective”  thresh¬ 
old  energies  leading  to  the  conclusion  that  the  “lucky  elec¬ 
tron"  picture  of  impact  ionization  is  inadequate  in  these 
important  semiconductors. 
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